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Abstract 

Quantum transitions between the Mott insulator and metals by controlling fill- 
ing in two-dimensional square lattice are characterized by a large dynamical 
exponent z = A where the origin of unusual metallic properties near the Mott 
insulator are ascribed to the proximity of the transition. The scaling near the 
transition indicates the formation of flat dispersion area due to singular mo- 
mentum dependence of the single-particle renormalization. The flat dispersion 
controls critical properties of the Mott transition. An instability of the flat dis- 
persion to the d-wave superconducting order is discussed. We also discuss a case 
of the Mott transition for a model of Mn perovskite compounds with orbital 
degeneracy where orbital correlation length shows critical divergence toward the 
metal-insulator transition. 

1 Introduction 

Metal-insulator transitions (MIT) driven by strong correlation effects are called 
Mott transitions and are a subject of recent intensive studies in two- and three- 
dimensional systems For one-dimensional (ID) systems, the Mott transitions 
are relatively well understood where rigorous results can be obtained for some 
nontrivial cases. In infinite-dimensional systems where the dynamical mean field 
theory becomes exact, the MIT can also be studied in a controlled way. How- 
ever, in ID systems, the Fermi level (and the Fermi surface) is represented by 
only two points while in the infinite-dimensional systems, the single-particle self- 
energy is site-diagonal and has no wavenumber dependence. Therefore, in both 
cases, charge excitations near the Fermi level are not allowed to have momen- 
tum dependence along the Fermi surface. On the contrary, as we see below, the 
momentum dependence of the single-particle renormalization may have singular 
dependence along the Fermi surface in finite-dimensional systems such as in two 
and three dimensions and therefore would show quite different features from one- 
and infinite- dimensional systems. 

In two or higher dimensional systems, there exists no evidence to exclude 
the Fermi liquid phase as the only one stable fixed point of ordinary metals. 
However, this is not necessarily an isotropic metal. When charge excitations in 
a part of the Fermi surface have slow dynamics, they are liable to couple strongly 
to other degrees of freedom ( such as spin fluctuations) thus yielding even slower 
dynamics. This synergetics may selfconsistently generate singular momentum 
dependence of the single-particle renormalization, where at a particular region 
(or points) on the Fermi surface, the Fermi liquid description breaks down. This 



1 



particular region may be characterized by flat dispersion with mass enhancement 
due to the real part of the selfenergy as well as by strong damping due to the 
imaginary part. 

Metallic phase near the Mott insulator is known to show various unusual 
properties in terms of standard metals. However, in any experimental systems 
clear phase boundary is not observed between the paramagnetic metallic phase 
near the Mott insulator and more or less standard metals observed far from 
the correlated insulator. This strongly suggests that the unusual properties are 
not due to the appearance of a new phase with adiabatic discontinuity from the 
weakly correlated metals, but due to a proximity of the MIT. Therefore, it is im- 
portant to understand critical properties of the MIT. Metal-insulator transitions 
are typical examples of quantum phase transitions which takes place at zero tem- 
perature by changing a parameter to control quantum fluctuations rather than 
thermal fluctuations. In this paper, we discuss that the unusual properties can 
be understood from novelty of the quantum criticality at the transition to the 
Mott insulator [|[ ^, In quantum phase transitions, dynamical fluctuations 
play roles even in high spatial dimensions. These are correctly treated in the 
"dynamical mean field theory" which provides exact result if spatial dimension- 
ality is infinite and spatial fluctuations are suppressed In real materials, 
however, quantum fluctuations with strong interaction effects are most conspic- 
uous in systems with low-dimensional anisotropy. In this case, fluctuations near 
quantum phase transition points are enhanced and appear as combined effects 
of strong spatial and dynamical quantum fluctuations. 

The Drude weight deflned from the coefficient of the (5-function in the frequency- 
dependent conductivity at T = as 

a{uj) — DS{uj) + crrcg(w) (1) 

and the charge compressibility k defined below are the two important and rel- 
evant quantities to describe the Mott transition. The Mott insulator has two 
basic properties: one is insulating property and the other, incompressibility. The 
Drude weight and the charge compressibility are both zero in the Mott insulat- 
ing phase while both nonzero in metallic phases. Therefore, they are totally 
nonanalytic at the transition point. The nonanalyticities in these two quantities 
are also common in the transition to the band insulator, while they do not have 
singularities in case of the Anderson localization transition. 

We here discuss in more detail the important difference between the Anderson 
transition and the MIT to the Mott insulator. In the Anderson transition, it is 
well known that the singularity at the transition point appears in the diffusion 
constant or the relaxation time of the carrier r (inverse of the imaginary part of 
the single-particle self-energy) but neither in the carrier effective mass to*, nor 
in the density of states N{uj) etc. Therefore the DC conductivity oc nr/m* is 
a good probe to determine the character of the transition because it is directly 
proportional to the relaxation time. In contrast, in the transition to the Mott 
insulator, the singularity may not appear in the relaxation time of the carrier 
but appears in the ratio of the carrier density n to the effective mass to*. The 
DC conductivity in the ideal case (without disorder at T = 0) is always infinite 
in the metallic phase and cannot be a good probe in the ideal condition (it can 
be a probe to some extent if the disorder or thermal effects play some role) while 
the Drude weight D oc n/m* and the charge compressibility k discussed here 
become relevant quantities with singularities. 

We note that the quantities which explicitly depend on the relaxation time 
arising from thermal or disorder effects have to be treated carefully. As we see 
below it is possible that the MIT is controlled by a particular part of flat disper- 
sion on the Fermi surface. However, the relaxation-time-dependent quantities 
such as DC transport properties can be "contaminated" by the relaxation time 
and the real criticality could be obscured because the contribution from the part 
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of flat dispersion could be suppressed due to shorter relaxation time as compared 
to the other part. 

We also note the difference between the Mott transition and magnetic quan- 
tum phase transitions which take place within metals. In the latter case, self- 
consistent renormalization theory was developed ^ where, at T = 0, a Gaus- 
sian treatment from the paramagnetic phase and the Hartree-Fock-RPA de- 
scription are justified above the upper critical dimension. There the anomalous 
character of metals is mainly ascribed again to the anomaly in the relaxation 
time although a weak effect on the renormalization factor is expected. 

In §2, we review results on novel universality class of two-dimensional MIT 
driven by correlation effects. The scaling theory with the universality class char- 
acterized by the dynamical exponent z = 4 offers a unified description for the 
unusual properties of metal near the Mott insulator. Numerical results on the 
charge compressibility, spin correlations, dynamical conductivity in the metallic 
phase as well as the localization length in the insulator are consistently un- 
derstood from this quantum criticality. The strong momentum dependence of 
the renormalization around (tt, 0) spot is suggested as the basis of the scaling 
description. Superconductor-insulator transitions have been studied motivated 
from instability of anomalous metallic states associated with this novel univer- 
sality class of the MIT. The results are summarized in §3. MITs with orbital 
degeneracy is the subject of §4. 



2 Metal-Insulator Transition in Two Dimensions 
2.1 Numerical Results 

The Mott insulator and metallic states near the MIT have been intensively 
studied subject by numerical approach. Here, we see that the 2D single-band 
Hubbard and t-J models show consistency with the scaling theory (|, It 
is discussed later that the single-particle excitation of the 2D Hubbard model 
around (±7r, 0) and (0, ±7r) in the momentum space is renormalized to form 
a flat quartic dispersion. This singular renormalization effect determines the 
universality class of the transition between the Mott insulator and metals and 
constitutes a basis for the scaling theory. This circumstance is a unique feature 
of the Mott transition in finite-dimensional systems in contrast with one and 
infinite dimensions. 

Before the clarification of the full momentum dependence, the singular mo- 
mentum dependence was suggested from the "momentum integrated" quantities. 
The first signature was observed in the charge compressibility. The charge com- 
pressibility K or charge susceptibility Xc are defined as Xc = n^n = dn/d^ where 
n is the electron density and fi is the chemical potential. The charge suscep- 
tibility vanishes in the Mott insulating phase because of the incompressibility. 
However, quantum Monte Carlo results on 2D systems show that doped systems 
become more and more compressible near the MIT in the metallic side |l^ . 
The scaling plot of Xc shows singular dependence on the doping concentration d 
in the form 

Xc cc l/\S\P (2) 

with p = 1 for (5 7^ 1^, |o[ |ll[. Similar is also observed in the 2D t-J model 
p^ . When the charge susceptibility is singularly divergent as a power of 6 
for S +0 in 2D, the single-particle description of low-energy excitations has to 
show the divergence of the effective mass of relevant particle. This is in contrast 
with the usual MIT between metals and the band insulator, where the number 
of carriers vanishes. 

The second signature of unusual renormalization was observed in the mea- 
surement of the localization length in the insulating phase [ p^ . A large numerical 
advantage to observe the transition from the insulator side is that we can avoid 
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the negative sign problem known in the quantum Monte Carlo method. The 
insulator undergoes a transition to a metal when the chemical potential fj, ap- 
proaches the critical point /ic from the region of the charge gap. In the insulating 
phase, the single-particle Green function defined as 

g{r,T) = {Tc{r,T)c^Q,0)) (3) 

provides the localization length defined by 

g{r,oj = n)= g{r,T,fi)dT^e-''/^'. (4) 

may be regarded as the localization length of the wavefunction of virtually 
created state at the chemical potential inside the gap. The scaling of ^/ near 
fic calculated with the above method by Assaad and Imada for the two- 
dimenional Hubbard model shows 

6 ~ l^J- - ^J-cl'" 

V = 0.26 ±0.05. (5) 

The obtained v is consistent with v — 1/4. 

The third signature is seen in the Drude weight. The Drude weight was 
calculated for the t-J model by the exact diagonalizationJTBI . The doping con- 
centration dependence in Fig. |l| shows that the total kinetic energy {~F) is 
proportional to 5 while the Drude weight D vanishes much faster consistently 
with D fx S^. It should be noted that this is in sharp contrast with a naive ex- 
pectation that the carrier density is proportional to S, where the Drude weight 
in the Drude theory is proportional to n/m* and should be linear in the carrier 
density n divided by the carrier effective mass m* . As we discuss in §2.3, the 
frequency dependence of the conductivity in numerical results also suggests that 
the most part of the conductivity weight in the sum rule is exhausted in a rather 
universal line shape of incoherent conductivity crreg(w) and is consistent with a 
quick collapse of the Drude part with decreasing doping. 

The above signatures observed in the charge compressibility, localization 
length and dynamical conductivity are hard to understand when isotropic renor- 
malization of charge excitations is assumed along the large Fermi surface which 
satisfies the Luttinger theorem. 

2.2 Scaling Theory of Mott Transition 

The numerical results on the charge compressibility and spin correlations in 2D 
have inspired extensive studies to clarify the nature of the Mott transition in 
low-dimensional systems H, ^ . Because the MIT is controlled by quantum 
fluctuations, we may take the control parameters such as the electron chemical 
potential /i or the bandwidth t. Using this control parameter, the distance from 
the critical point is measured by A. The control parameter can either be the 
chemical potential to control the filling or the bandwidth (or the interaction). 
The scaling theory assumes the existence of single characteristic length scale ^ 
which diverges as |A| and a single characteristic frequency scale H. which 
vanishes as |A| ^0. The correlation length ^ is assumed to follow ^ 1^1^'' as 
A ^ 0, which defines the correlation length exponent i'. The frequency scale 
is determined from the quantum dynamics of the system independently of the 
length scale. The dynamical exponent z determines how vanishes in relation 
to 1/^ as ri ~ ^ jAj^". The hyperscaling asserts the homogeneity in the 
singular part of the free energy density fg in d + I dimensional path integral in 
terms of arbitrary length-scale transformation parameter b: 

/,(A) ~ b-'^'^+'^f,ib^/''A) ~ A'''-'^+'\ (6) 
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When the inverse temperature (3 and the hnear dimension of the system size L 
are both large but finite, finite-size scahng function is expected to hold as 

/,(A)^A''(''+^)^(e/L,e7/3) (7) 

in the combination of non-dimensional arguments. 

Combining expressions in the path integral formalism for the Drude weight 
and the compressibility, respectively, with the finite-size scaling form derived 
from the hyperscaling, one obtains useful scaling forms for physical quantities 
0, H . The scaling of D is obtained: 



with ( — v{d + z — 2), while 
and 



D oc A'^ (8) 

■ cx A-"+''" (9) 
KocA"" (10) 



are obtained with a = v{z — d), where 5 denotes the carrier concentration, 
namely, the density measured from the Mott insulating phase. 

In the insulating phase, the charge excitation gap Eg and the localization 
length defined in should be determined from Eg oc lAj^" and 

« |Ar^ (11) 

where we have used the fact that v and z are identical in both sides of the 
transition point. In the metallic side, the scaling of the specific heat C = 
f3^d^{liiZ)/d(3^ is given from (g): 

C = A'-^'+^^T-^TiC/L = 0, TT). (12) 

At the critical point, the specific heat is given from the A independent term as 
C cx T"^/^, If C* = is satisfied at low temperatures in the metallic phase, the 
coefficient 7 follows 

7cxA''(''-^^ (13) 

Another interesting quantity is the coherence temperature Tp below which the 
electron motion becomes quantum mechanical and degenerate. It is clear that 
Tp has to approach zero as A —> in a continuous transition. In case of the 
Fermi liquid, Tp is nothing but the Fermi temperature. The existence of single 
characteristic energy scale with singularity at A = leads to the scaling of Tp 
in the form 

Tp (X A"'. (14) 

When the control parameter A is the chemical potential fi measured from 
the critical point /ic, it represents the filling control MIT. In this case, it is 
possible to derive a useful scaling relation. Because the doping concentration is 
5 = = — scales as 

From the comparison ofEq.(|l|) and Eq.(|), j^z = 1 is derived. The characteristic 
length scale is then 

e^,5-i/^ (16) 

which is the length scale of the "mean hole distance" . The above scaling de- 
scription is known to be valid for the MIT to the band insulator, all the MIT in 
ID [|r^ and the Anderson transitions ]l8| |. 



5 



2.3 Universality Class 

When we assume the above scahng theory with Eg. ([l0|) , the quantum Monte 
Carlo resuh (||) imphes that there exists a new universality class z — 1/ v — 4^ 
[H, ^. As we see above in (^) and (p]), the locahzation length also suggests 
that V ^ 1/A. These are independent estimates of the exponent and a check for 
the consistency of the scaling theory. On the contrary, these two results can be 
explained neither by the Hartree-Fock approximation nor by the d = oo results. 
Indeed, neither the Hartree-Fock nor the d = oo results satisfy the hyperscaling 
assumption (^. For example, the Hartree-Fock approximation predicts v = 1/2. 

The large dynamical exponent z leads to unusual suppression of coherence 
in the metallic phase in various aspects. The scaling theory predicts that the 
Drude weight scales as 13 oc 5^'^ ~ . As is well known, 13 oc (5 is satisfied in 
the usual MIT, which is consistent with the scaling theory when z = 2. When 
the universality class is characterized by z = 1/z/ > 2, 13 is suppressed stronger 
than 5-linear dependence at small 5. When z = A is satisfied in 2D, D o^ 5^ is 
predicted. This is one of the indications for the unusual suppression of coherence 
in this universality class. Numerical results indeed showed consistency with 
Z3 cx (5^ as we discussed in §2.1 [ p^ . 

From the sum rule, the tj-integrated conductivity gives the averaged kinetic 
energy 

/>oo 

~ {F)^ I a{uj)du. (17) 

In the strong coupling limit as in the t-J model, the total kinetic energy {F) is 
expected to be proportional to 5 as observed in Fig. la. Then the Drude part 
D oc (5^"' ~ becomes negligibly small at small 5 in the total weight —(F) oc (5 if 
z > 2. This means most of the total weight —{F) is exhausted in the incoherent 
part (Trcg(u;) in (|l|). In the scaling theory, the form of (7rcg('^) is not specified and 
indeed it may depend on details of systems including the interband transition. 
However, for the intraband contribution, a{u}) is expected to follow <j{uj) 
(1 — e^^'^)/ui at very high temperatures T larger than the bare bandwidth Ib- 
In fact the conductivity is given by 

1 - e'^P'^ 

a{oj) = C{uj) (18) 

with the current correlation function C{ljj) defined by 

C{u) ^ r dt e^-'{mm) ^ (19) 

where if the carrier dynamics is incoherent, C{uj) is expected to have a broad 
featureless structure due to rapid decay of current correlation in time t l/ts- 
When the temperature is lowered, this incoherent part ~ (1 — e~^^)/uj is in 
general transferred to the Drude part ~ below the coherence temperature. 

In the usual band-insulator-metal transition, this transfer can take place for the 
dominant part of the weight because — {F) and D both may be proportional to 
5. Therefore, the incoherent part is totally reconstructed to the Drude part with 
decreasing temperature below Tp. However, for z > 2, this transfer takes place 
only in a tiny part of the total weight because the relative weight of the Drude 
part to the whole weight is ij^^"^)/'' which vanishes as (5 — > for z > 2. Even 
at T = 0, the major part of the conductivity weight must be exhausted in the 
incoherent part. Near the transition point, since we have no reason to have a 
dramatic change of the form of the incoherent part o'iog('i^) at any temperatures 
from the scaling point of view, the optical conductivity more or less follows the 
form 

1 _ p~P^ 

CTreg(w) - C (20) 
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even for > and ts > T where we have neglected model dependent feature 
such as the interband transition. This "intraband" incoherent weight may be 
proportional to C oc in the strong coupling limit. Such dominance of the 
incoherent part with the scaling (EO) is consistent with the numerical results of 



2D systems 19 , which lends further support for 2; > 2 in 2D. On the contrary, the 
numerical results in ID by Stephan and Horsch |2^ support that the incoherent 
part is small at low temperatures. This may be due to z = 2 in ID because the 
majority of the weight seems to be exhausted in the Drude weight. 

Another indication for the suppression of coherence may be seen in the co- 
herence temperature Tp. The scaling of Tp is given by Tp oc 5^^'^. Because the 
standard MIT is characterized by TpQ cx (5^/'^, the relative suppression Tp/Tpo 
is proportional to S~3- _ When z = 4 as suggested by numerical results in 2D, 
we obtain Tp oc S'^ in 2D in contrast with TpQ oc S for z = 2. 

In the Mott insulating phase, entropy coming from the spin degrees of free- 
dom is released essentially to zero when the antiferromagnetic order exists. The 
entropy is also released when the spin excitation has a gap as in the spin gap 
phase. Even without such clear phase change, the growth of short-ranged cor- 
relation progressively releases the entropy with decreasing temperature. When 
carriers are doped, an additional entropy due to the charge degrees of freedom 
is introduced in proportion to 5. This additional entropy is assigned not solely 
to the charge degrees of freedom but also to the spin entropy through their 
mutual coupling which destroys the antiferromagnetic long-range order. This 
additional entropy oc d has to be released below the coherence temperature Tp. 
Therefore, a natural consequence is that if T-linear specific heat characterizes 
the degenerate (coherent) temperature region, the coefficient 7 should be given 
by 7 ^ S/Tp ^ 6^^^. This is indeed the scaling law we obtained in ( [l3|) and 
(|l5|). From this heuristic argument, it turns out that the metallic phase near the 
Mott insulator is characterized by large residual entropy at low temperatures 
which is also related with the suppression of the antiferromagnetic correlation 
as compared to the insulating phase. It may also be said that the anomalous 
suppression of the charge coherence at small S is caused by short-ranged anti- 
ferromagnetic correlations which scatters carriers incoherently. 

The strong coupling of spin and charge seems to yield the same scaling even 
for the spin correlation. The numerical result suggests that the antiferromagnetic 
transition takes place more or less simultaneously with the Mott transition ||^. 
The equal time spin structure factor S{Q) at its maximum shows the scaling 
S{Q) oc 1/6 from which the magnetic correlation length ^.afm is scaled by 
^-1/2^ The incommensurate wavevector Q approaches (7r,7r) as 5 — > 0. This is 
the same scaling as (1™ and implies that the criticality of the antiferromagnetic 



transition is also involved in the same universality. 



2.4 Flat Dispersion and Strong Momentum Dependence 

To clarify the origin of the large dynamical exponent, z = 4, it is necessary to 
study explicit momentum dependence of the charge excitations. The realization 
of the hyperscaling must be the consequence of the fact that some small number 
of singular points in the momentum space appear and the MIT is controlled by 
the charge excitations around those points. If the charge excitations near the 
whole d — 1 dimensional Fermi surface would contribute equally and isotropi- 
cally, it would not satisfy the hyperscaling at all because it would introduce a 
length scale of the Fermi wavenumber which does not vanish at the MIT point. 
This additional length scale would invalidte the assumption (^). Furthermore, 
the dynamical exponent z = 4 directly has to lead to the appearance of a fc^ 
dispersion around those singular points because by definition z is nothing but 
connecting characteristic wavenumber and energy scales. If a flat dispersion 
such as fc* appears at some points on the Fermi surface, the holes will be doped 
predominantly to this region because of this flatness and this flat region will 
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determine the character of the MIT. 

To understand whether the above reasoning works, the single-particle spec- 
tral function A(k,uj) was calculated and analyzed in detail. A{k,uj) of a doped 
single hole to the Mott insulator phase of the Hubbard model has peak struc- 
tures which disperses in the momentum space. Figure |^ shows this dispersion 
around (tt, 0) and clearly indicates that the flat dispersion does appear around 
(tt, 0) and (0, tt), while other part of the Fermi level does not show such flat- 
ness. The numerical result further indicates that the dispersion around (tt, 0) 
and (0, tt) is well fit by fc^ dispersion. These points have originally the van-Hove 
singularity in the simple Hubbard model. However, the observed k'^ dispersion is 
far beyond the expectation from the van-Hove singularity and the consequence 
of strong correlation effects. We further note that this z = 4 universality class 
also holds for the case with finite next-nearest-neighbor hopping t' where the 



van-Hove singularity does not meet the Fermi level at the MIT point |10 . This 



gives us the microscopic basis to understand the anomalous feature of the MIT. 

The dispersion of a doped hole in the Mott insulating state does not neces- 
sarily determine the criticality of the transition because the rigid coherent band 
picture is not always guaranteed in the process of further doping. However, the 
existence of the flat dispersion is commonly observed at low doping concentra- 
tion in angle resolved photoemission (ARPES) data of high-Tc cuprates pl|, |2^. 
The presence of this flat dispersion is also consistent with numerical results of 
the Hubbard models at finite doping |Q . More remarkable is the following: The 
ARPES data for the insulating phase of the cuprates appears to show a deeper 
level around (tt, 0) than the level around {tt/2, 7r/2) which is relatively closer to 
the Fermi level as is consistent with the case with the next nearest neighbor 
transfer t' [Q. If the rigid band picture would be correct, this implies that 
the state around (7r/2,7r/2) would be first doped with holes rather than (7r,0) 
region. However, in the underdoped region, the data indicate that the state 
around (tt, 0) seems to quickly emerge from the incoherent tail of A{k, uj) at the 
Fermi level and forms a flat dispersion, which indicates the breakdown of the 
rigid band picture and a universal dominance of the flat dispersion excitations 
upon doping. This is also consistent with the numerical compressibility data 
with t' discussed above. The fiat dispersion could be obscured by its incoherent 
nature and hard to observe incoherent tails in A{k,uj) by the strong damping. 

The dominance of the flat dispersion around (tt, 0) on low energy excitations 
should be taken with care because it has some degeneracy with the excitations 
around (7r/2,7r/2). It may have an effect on the DC transport properties, al- 
though the DC properties are not good probes to understand the Mott transi- 
tion as we discussed in the introduction. The ARPES data of the underdoped 
cuprates show a gradual formation of fragmentary "Fermi surface" first around 
(7r/2,7r/2) in the normal state. This may have two reasons. One is that the 
carrier relaxation time around (7r/2,7r/2) is much longer and the other is that 
the region around (tt, 0) is under the influence of the formation of pseudogap. 
The low-energy excitation around (vr, 0) is dominated by the paired singlet for- 
mation due to the instability of the flat dispersion as we will discuss in the next 
section. The paired singlet formed in the pseudogap region seems to be still 
incoherent above Tc |2^, Although the MIT is governed by the excitations 
around (tt, 0), the DC transport may have substantial contribution from the part 
around {tt/2, 7r/2). This is speculated from the fragmentary Fermi surface as well 
as from a rather large and sensitive increase of the residual resistivity upon Zn 
doping around 1% in the cuprates while relatively insensitive effect of Zn on the 
pseudogap behavior ||2^, ^ . This implies that the Zn doping sensitively induces 
the localization of carriers around (7r/2,7r/2), while gapped singlet excitations 
formed from (tt, 0) excitations is not seriously influenced by such tiny amount 
of Zn. The DC charge transport is strongly affected from (7r/2,7r/2) excitations 
while other quantities including spin excitations are dominated by the driving 
excitations of MIT near (7r,0), which makes the physical properties as if spin 
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and charge degrees of freedom were "separated" . However, the seeming separa- 
tion has to be understood after considering the strong momentum dependence 
carefully. 

3 Superconductor-Insulator Transition 

In the Mott insulator, because the single-particle process is suppressed due to 
the charge gap, we have to consider the two-particle process expressed by the 
superexchange interaction. This has led to the Heisenberg model for the effec- 
tive Hamiltonian of the Mott insulating phase. The Hubbard model would be 
sufhcient even for the Mott insulating phase if its low energy excitations could be 
precisely considered. However, the Heisenberg model offers a better and easier 
way to extract physics of spin excitations and properties at low temperatures in 
the strong coupling regime. Similarly to this circumstance in the insulating case, 
even in the metallic region, such suppressions of coherence discussed above for 
the single-particle process make it useful and helpful to consider the two-particle 
process explicitly. This is particularly true in the flat-dispcrsion part of the mo- 
mentum space. In contrast with the insulating phase, the two-particle processes 
in metals contain an effective pair hopping term and a term described by 

Uw = -tw^^{c\^c^+sa + (21) 

i 5a 

becomes relevant ^ Then the Hamiltonian containing the term TLw 
added to the Hubbard model, called the t-U-W model is expected to be a useful 
effective Hamiltonian to understand low temperature properties near the Mott 
insulating phase. It was argued that Tiw plays the same role as the term pro- 
portional to /U in the strong coupling expansion of the Hubbard model which 



contains the superexchange term as well as the so-called three-site term 30 . It 
is not clear at the moment whether the bare Hubbard model implicitly contains 
sufficiently large effective coupling of Tiw to make it relevant at low energy scale 
or some additional elements are needed. However, by including a small amount 
of the TLw term explicitly, it allows to study how the Hubbard hamiltonian be- 
comes unstable to the two-particle process and what type of symmetry breaking 
or ordering are expected as its consequence. In fact the t-U-W model is the first 
hamiltonian which has made possible to study the d-wave superconductor-Mott 
insulator transition under a controlled numerical treatment. 

We note that, to treat the relevance of the Tiw term under a proper condition, 
it would be necessary to first extract precisely the strong momentum dependence 
of the single-particle renormalization in the Hubbard model. Without clarifying 
the formation of the flat dispersion caused by many-body effects described in the 
previous section, the real role and relevance of the two-particle term cannot be 
understood well. Mean field analyses of the Hubbard as well as the t-J models 
in the literature have not properly treated this important point by neglecting 
the singular momentum dependence of the renormalization. 

When the two-particle process becomes relevant, the coherence is not sup- 
pressed because the universality class of the two-particle transfer is given by 
z = \/v = 2. The universality class z = 2 is numerically observed in the expo- 
nent of the localization length in the insulating side of the t-U-W model ||3^ and 
the doping concentration dependence of the Drude weight of the t-J-W model 
[ p^ . Because of small z, it always becomes more relevant than the single-particle 
transfer for small 5. The alteration of the universality class from = 4 for the 
MIT to z = 2 for the superconductor-insulator transition shows an instability of 
metals near the Mott insulator to the superconducting pairing or in more gen- 
eral an instability of z = 4 universality phase to some type of symmetry-broken 
state. 
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Below, we summarize how the Mott insulator to superconductor transition 
takes place and how the d-wave superconducting state appears in the 2D t-U- 
W model. This has been studied through large-scale quantum Monte Carlo 
calculation ||2^, |3^. The t-U-W model appears to show a Mott insulator to 
superconductor transition even when the filling is fixed at half filling. With the 
increase of W, the antiferromagnetic long-range order seen in the pure Hubbard 
model continuously decreases and at the critical amplitude oi W = Wc, the 
order is destroyed. At the same W, the insulating phase appears to undergo a 
transition to a superconductor with d^2_y2 symmetry. 

A remarkable property in the superconducting phase at half filling is that 
the antiferromagnetic correlation is extremely compatible. The equal-time anti- 
ferromagnetic correlation decays very slowly with a power law at long distance 
as ^ l/r" with a ~ 1 or slightly larger. Reflecting this compatibility, the real 
part of the staggered susceptibility xil = {^^i'^):^ = 0) appears to diverge with 
lowering temperature to the limit T — > 0. A clearer understanding is obtained 
in the dynamical structure factor S{q,uj) (or equivalently in the imaginary part 
of the dynamical susceptibility !mx{q,uj))- From this quantity, it turns out that 
the antiferromagnetic correlation is dynamical and has a strong peak at a fi- 
nite frequency with presumable divergent weight in the thermodynamic limit. 
This peak frequency is comparable to the pairing gap amplitude. Because of 
the appearance of the superconducting phase, the antiferromagnetic correlation 
is suppressed below the pairing energy scale. However very slow power-law de- 
cay of equal-time antiferromagnetic correlation implies that S{q = (7r,7r),w) has 
divergent weight somewhere at finite frequency. To reconcile the antiferromag- 
netic and superconducting correlations, the frequency of the antiferromagnetic 
correlation is pushed out from low-frequency to the frequency above the pairing 
gap. When the temperature is lowered at Vl^ > Wc from high temperatures, 
S{q = (tt, tt), lo) first shows growth of a broad peak around lu — 0. This appears 
to be an immatured state which does not differentiate the antiferromagnetic and 
singlet correlations. This is reminiscent of the S0(5) scenario with approximate 
symmetry of superconducting and antiferromagnetic state , although the ap- 
proximate symmetry in this model has not been examined in detail yet. With 
further lowering of the temperature, the peak position shifts to a finite frequency 
and is sharpened. In this temperature range pseudogap 

is formed where lmx{q = (tt, 7r),a;) decreases at low lu. At lower tempera- 
tures, the peak at finite frequency becomes even sharper. 

4 Mott Transition with Orbital Degeneracy 

Orbital degeneracy plays important roles in addition to the spin degeneracy in 
a wide class of materials at the MIT [Q . It provides another origin of residual 
entropy near the transition point and may offer richer structure of the transition. 
Effects of orbital degeneracy was recently examined by projector quantum Monte 
Carlo method at zero temperature for a model of Mn perovskite compounds with 
degeneracy of two eg orbitals, dj.i_y2, and d^z'^-r'^, under the condition of com- 
plete spin polarization and two-dimensional configuration p3[ . In Mn perovskite 
compounds, the Mott insulator is in the configuration where one of two Cg 
orbitals are occupied. Through the MIT, the spin configuration is always com- 
pletely ferromagnetic within a plane mainly due to a strong Hund's rule coupling 
to high-spin t2g electrons. However, the charge transport is strongly incoherent 
with very small Drude weight in the metallic region js^ . This suggests a crucial 
role of orbital degeneracy (and presumably also the lattice degrees of freedom 
coupled to it) for the incoherent charge dynamics. The spin polarized but or- 
bitally degenerate model may be a minimal model to capture this circumstance. 
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The Hamiltonian of this model is defined as 

+U ^ {ni^ -^){n^^' - fJ.^n^„, (22) 

where the creation (annihilation) of the single-band electron at site i with orbital 
V is denoted by cl^ici,^) with rii^ being the number operator riii, = cl^Civ An 
important difference from the ordinary Hubbard model is that the transfer has 
an anisotropy with dependence on the — and — orbitals, where we 
label the orbital — and 3z^ — as 1 and 2, respectively below. In 2D 
configuration, the nearest-neighbor transfer is scaled by a single parameter t as 
til — jt, t22 = jt, ti2 ~ <2i = i-^^j where, in ±, +(— ) is for the transfer 
to the y{x) direction. In the absence of the diagonalized two bands have 
dispersions both with the bandwidth 4t, where one dispersion is obtained by a 
parallel shift of 2t from the other. 

The present numerical result is summarized as follows. When U is increased, 
the Mott gap Ac opens at "half filling", namely at (n) — 1. However, the 
opening of the gap is substantially slower than the case of the ordinary single- 
band Hubbard model (with spin). For example, Ac ~ 0.1 at U/t = 4 is compared 
with Ac ~ 0.66 for the ordinary Hubbard model at U/t = 4. At U/t — 3, Ac is 
not distinguished from zero within the error bar. Although the opening of the 
Mott gap is slow within this model, the insulating state at (n) = 1 is stabilized 
if we introduce the Jahn Teller distortion. A realistic Jahn- Teller coupling for 
the Mn compounds has stabilized the Jahn- Teller distorted insulating state with 
staggered order of — r^ and 3?/^ —r^ orbitals. Another important observation 
is that realistic amplitudes of the Jahn- Teller coupling in the absence of U is far 
insufficient in stabilizing the Jahn- Teller distortion with a realistic stabilization 
energy of several hundred K. The Mott insulating state in the experimental 
situation may well be resulted from synergy of U and the Jahn- Teller coupling. In 
the model ( p2|) with 2D configuration, the long-ranged orbital order increases at 
(n) = 1 with increasing U . For example the staggered orbital polarization at [/ = 
4t is around 0.3. However, the orbital order appears to be destroyed immediately 
upon doping. In the metallic region n ^ 1, the short-ranged correlation of the 
staggered orbital order is critically enhanced as (n) — > 1 similarly to the ordinary 
2D single-band Hubbard model. For large J7, there seems to exist a critical region 
near (5 = 1 — n = where T(ff) oc 1/5 as the same as the spin structure factor 
of the ordinary Hubbard model. Here T{Q) is the equal-time structure factor 
at (tt, it) for the orbital correlation. This critical enhancement of the orbital 
correlation strongly suggests that the orbital correlation in this Mn model may 
play a similar role to spins in the Hubbard model in realizing highly incoherent 
charge dynamics. 
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Figure 1: (a) Drude weight D and total weight {—F) in the optical conductivity 
of the 2D t-J model with 4x4 sites at zero temperature, (b) The ratio of the 
two weights D/{-F). 
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Figure 2: The dispersion of the Hubbard model along the kx axis, from (0,0) 
through (tt, 0) till (27r, 0) obtained from the peak position of the single-particle 
spectral weight A{k, oj) for the 16 by 16 lattice at zero temperature. Black circles 
are the Monte carlo data at U/t — 4. 
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